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Within the framework of non-standard (Weyl) representations of the canonical commutation re- 
lations, we investigate the polymer quantization of the Taub cosmological model. The Taub model 
is analyzed within the Arnowitt-Deser-Misner reduction of its dynamics, by which a time variable 
arises. While the energy variable and its conjugate momentum are treated as ordinary Heisenberg 
operators, the anisotropy variable and its conjugate momentum are represented by the polymer tech- 
nique. The model is analyzed at both classical and quantum level. As a result, classical trajectories 
- - - , flatten with respect to the potential wall, and the cosmological singularity is not probabilistically 

9^ ' removed. In fact, the dynamics of the wave packets is characterized by an interference phenomenon, 

' which, however, is not able to stop the evolution towards the classical singularity. 

^ , PACS numbers: 98.80.Qc;11.10.Nx 

> 
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■ I. INTRODUCTION 

The necessity for a quantum theory of gravity arises from fundamental considerations, and, in particular, from the 
O space-time singularity problem. In fact, the classical theory of gravity implies the well known singularity theorems, 
^ '' among which the cosmological one The canonical quantization of gravity, which exhibits a host of difficulties 
5—1 , both at technical and interpretative levels, is based on the Heisenberg representation of the Weyl algebra On the 
, other hand, the background-independent formulation of canonical quantum gravity based on Yang-Mills formalism 

has recently appeared ^3] • Anyhow, smearing such variables in the holonmy-flux representation is an important step 
PsJ ■ towards canonical quantum gravity [3| ■ The scenario induced by such an algebra is illustrated to be equivalent to the 
J> ] so-called polymer representation of quantum mechanics [ESB]' ^ soon as a mechanical system is taken into account. 
,92 ■ This work is aimed at investigating the quantization of the Taub model in the polymer representation of quantum 
mechanics. 

The Taub Universe arises as a particular case of the Bianchi IX model, i.e. the most general scheme allowed by 
the homogeneity constraint @. In the Bianchi IX model, the Universe dynamics towards the classical singularity is 
summarized by the chaotic motion of a particle. More precisely, this particle bounces an infinite number of times 
against the potential walls of a triangular domain, on a two-dimensional plane. The two-dimensional plane describes 
^ the configuration space of the particle (Universe) dynamics. The Taub model consists in restricting the dynamics to 
• • that of a one-dimensional particle bouncing against a wall, when only one degree of freedom is taken into account. 
. ^ The relevance of the Taub universe in quantum cosmology is due to the fact that it is a necessary step towards the 
, more general Bianchi IX model. The advantage of this model is that it is a generalization of other isotropic models. 
^ • In particular, it has been used to test the validity of the minisuperspace scheme @ and to explore the application of 
5^ ] the extrinsic cosmological time [l^. Furthermore, the Taub model has also been investigated within the framework of 
a generalized uncertainty principle in where the cosmological singularity has been shown to be probabilistically 
removed. 

The polymer representation of quantum mechanics is based on a non-standard representation of the canonical com- 
mutation relations ff>]. In particular, in a two-dimensional phase space, it is possible to choose a discretized operator, 
whose conjugate variable cannot be promoted as an operator directly. From a physical point of view, this scheme 
can be interpreted as the quantum- mechanical framework for the introduction of a cutoff. Its continuum limit, which 
corresponds to the removal of the cutoff, has to be understood as the equivalence of microscopically-modified theories 
at different scales [7]. 

This approach is relevant in treating the quantum-mechanical properties of a background-independent canonical 
quantization of gravity. In fact, the holonomy-flux algebra used in Loop Quantum Gravity reduces to a polymer-like 
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algebra, when a system with a finite number of degrees of freedom is taken into account [5| . From a quantum-field 
theoretical point of view, this is substantially equivalent to introducing a lattice structure on the space [l^- Loop 
Quantum Cosmology can be regarded as the implementation of this quantization technique in the minisuperspace 
dynamics fT3|. 

The Taub model is approached in the scheme of an Arnowitt-Deser-Misner (ADM) reduction of the dynamics in the 
Poincare plane. As a result, a time variable naturally emerges, and the Universe is described by an anisotropy-like 
variable. The anisotropy variable and its conjugate momentum are quantized within the framework of the polymer 
representation. More precisely, the former appears as discretized, while the latter cannot be implemented as an oper- 
ator in an appropriate Hilbert space directly, but only its exponentiated version exists. The analysis is performed at 
both classical and quantum levels. The modifications induced by the cutoff scale on ordinary trajectories are analyzed 
from a classical point of view. On the other hand, the quantum regime is explored in detail by the investigation of 
the evolution of the wave packets of the universe. 
Two main conclusions can be inferred. 

• An interference between the wave packets and the potential wall appears. Nevertheless, the classical cosmological 
singularity is not probabilistically removed. In fact, the wave function of the universe is not strictly localized 
away from it, and the wave packets fall into it following a classical trajectory. 

• The comparison between the polymer approach and the Generalized Uncertainty Principle (GUP) model il- 
lustrates that the corresponding interference phenomena are produced in a complementary way. This feature 
appears both at classical level, as it is immediately recognized analyzing the modifications of the equations of 
motion, and in the quantum regime, as the behavior of the wave packets is investigated. 

The paper is organized as follows. 

In Section II, we review the main features of the Taub cosmological model. 

In Section III, the polymer representation of quantum mechanics is developed at both kinematic and dynamical level. 
Furthermore, the continuum limit of this approach is discussed. 

The fourth Section is devoted to the application of the polymer paradigm to the Taub universe, at both classical and 
quantum level. 

Section V is aimed at constructing suitable wave packets and investigating their dynamics. In particular, they are 

analyzed in the WDW and the polymer representation. 

In Section VI, our results are discussed and compared with other models. 

Concluding remarks follow. 

Throughout the paper, we have adopted natural units, i.e. h — c — IGttG — 1. 

II. THE TAUB MODEL 

Homogeneity reduces the configuration space of General Relativity to three dimensions. The homogeneous cosmo- 
logical models the Bianchi Universes, are such that the symmetry group acts simply transitively^ on each spatial 
manifold. The Bianchi IX model, together with Bianchi VIII, is the most general one and its line element reads, in 
the Misner parametrization jl5i] . 

ds"^ ^ N^dt^ -e'^'' {e^-')^^uj'(S,uj^, (1) 

where TV — N{t) is the lapse function and the left invariant 1-forms — oj^^dx"' satisfy the Maurer-Cartan equation 
2du;* = ejf.uj^ Auj'^ . The variable a — a{t) describes the isotropic expansion of the Universe and 7^ = 7y (*) is a traceless 
symmetric matrix, 7^ = diag (7-1- -f y/Sj-, 74. — \/37-, , which determines the anisotropy changes via 7±. The 

classical singularity appears for a ^ — 00, since the determinant of the 3-metric is given hy h = dete'^^'^'^ = e^". 

a. Canonical Analysis. The Hamiltonian constraint for this model is obtained performing the usual Legendre 
transformation. As well-known (isl . [l6| . the dynamics of the Universe towards the singularity is described by the 
motion of a two-dimensional particle (the two physical degrees of freedom of the gravitational field) in a dynamically- 
closed domain. Such a domain depends on the time variable a in the Misner picture, while it is stationary in the 



^ Let G a Lie group, G is said to act simply transitively on the spatial manifold E if, for all p, g G S, there is a unique element g £ G such 
that g{p) = q. 



FIG. 1: The dynamical-allowed domain Tq{u,v) in the Poincare complex upper half- plane where the dynamics of the Universe 
is restricted, towards the classical singularity, by the potential. 



Misner-Chitre framework defined by the variables [17 1 



Q- = -e^^, 7+ = eWi^ - lcos6l, 7_ = e^^C^ - lsin6', (2) 

witli ^ S [l,c») and 6 € [0, 27r]. In fact, the dynamically-allowed domain becomes independent of r, which behaves 
like a time variable. In terms of these new variables, the Hamiltonian constraint rewrites 

H = -pl+plie-l) + ^^^Q. (3) 

b. ADM Reduction. Let us perform the ADM [3| reduction of the dynamics. This scheme relies on the idea to 
solve the classical constraint with respect to a given momentum, before implementing any quantization algorithm. 
This paradigm allows us to dynamically separate the six-dimensional phase space of the model. In particular, a time 
variable arises and an effective Hamiltonian, which will depend only on the physical degrees of freedom of the system 
(the anisotropy-like variables), comes out. We solve explicitly the constraint H — with respect to Pr, and thus we 
consider the variable t as the time coordinate for the dynamics (we adopt the time gauge f = 1), obtaining 



-Pr = ^jplie-1) + ^^- (4) 

The dynamics of such a system is equivalent to a billiard ball on a Lobatchevsky plane , as we can see by means 
of the Jacobi metric^ . It is possible to choose the so-called Poincare representation in the complex upper half-plane 
[20| by using new variables {u, v), defined as 

^^l±J^±^^ e = -t.n-^( f^^ty (5) 
^ ^/3v y^l + 2u + 2u^ + 2v^ J 

The dynamical-allowed domain Fg — Tq{u,v) is plotted in Fig. 1. It is worth noting that the three corners in the 
Misner picture are replaced by the points (0,0), (—1,0) and w — > oo in the {u — v) plane. In this scheme, the ADM 
"constraint" is simpler than the previous one and becomes 



- = Hadm = v\/pl +pI. (6) 

The Taub Universe corresponds to the Bianchi IX one in the particular case of 7_ =0 Q . The phase space of this 
model is four-dimensional and its dynamics is equivalent to the motion of a particle in a one-dimensional domain. 
Considering such a domain corresponds to taking only one of the three equivalent potential walls of the Bianchi IX 



^ This approach reduces the equations of motion of a generic system to a geodesic problem on a given manifold. 
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model. As we can see from ([2]) and (O, this particular case appears for 6* = ^ u = -1/2 = (w^ + 3/4)/V3w), and 
the ADM Hamiltonian ^ rewrites 

Hadm = vpv, (7) 

being v G [1/2, cxd), as shown in Fig. 1. The Hamiltonian above ([7]) can be further simplified defining a new variable 
X = hiv, and becomes 

Hadm =Px=P, (8) 

which will be the starting point of our analysis. Within this framework, the Taub model is therefore described by a 
two-dimensional system in which the variable r is considered as the time, while the variable x describes the single 
degree of freedom of the Universe, i.e. the shape change. It is worth stressing that the classical singularity now 
appears for r — > cxd. 



III. POLYMER QUANTUM MECHANICS 



The polymer representation of quantum mechanics consists in defining abstract kets, labeled by a real number, 
and then considering a suitable finite subset of them, whose Hilbert space is defined by the corresponding inner 
product . This procedure can be shown to be an inequivalent representation of the Weyl algebra wrt the ordinary 
Schroedinger one. This representation helps one gain insight onto some particular features of quantum mechanics, 
when an underlying discrete structure is somehow hypothesized. The request that the Hamiltonian associated to 
the system be of direct physical interpretation defines the polymer phase space, and the continuum limit can be 
recovered by the introduction of the concept of scale Q • 

c. Wave functions and operators. One can start by considering abstract kets |^ >, ^ G R, and a suitable subset 
defined by /^i € R, i = 1,2, ..N. These kets are assumed to be an orthonormal basis, i.e., < fj,\v >= 6fj,u, along which 
any state can be projected. This defines a Hilbert space Tipoi, on which two basic operators act, the symmetric 
"label" operator, e, such that e|/i >= fi\fi >, and a one-parameter family of unitary operators, s(A), such that 
s(A)|/i >= \fi + X >. Because all kets are orthonormal, s(A) is discontinuous, and cannot be obtained from any 
Hermitian operator by exponentiation. It is worth noting that this Hilbert space is not separable^. 

For the toy model of a 1-dimensional system, whose phase space is described by the variables p and q, the polymer 
representation techniques find interesting applications when one of the two variables is supposed to be discrete. This 
discreteness will affect both wave functions, obtained by projecting the physical state on the p or q basis (polarization), 
and the operators associated to the canonical variables, acting on them. 

For later purposes, we will discuss only the case of a discrete position variable q, and the corresponding momentum 
polarization. 

In this case, wave functions are given by V'm(p) —"^ pIa* >= e'^''. Accordingly, the "label" operator i is easily identified 
with q, i.e., q4>f^ = —idpij}^ = fJ-ipfj., while the "shift" operator does not exist, as discussed previously. 
It can be shown that the corresponding Hilbert space is Hpoi = £^(Rb, d^//), i.e. the set of square-integrable functions 
defined on the Bohr compactification of the real line Rs, with a Haar measure dfin- Since the kets |/i > are arbitrary 
but finite, the wave functions can be interpreted as a quasi-periodic function, with the inner product 

<'0mIV'a>= / d^lH■ip^L{p)'^P\{p) = }^rn -!- [ dp^/j^{p)'il;x{p) = Sf,.x. (9) 

d. Dynamics. The Hamiltonian operator H describing a quantum-mechanical system is usually a function of 
both coordinate an momentum, i.e. H = H{q,p) = + V{q) while, in the particular case of a discrete position 
variable in the momentum polarization, p cannot be implemented as an operator, so that some restrictions on the 
model have to be required. 

As a first step, a suitable approximation for the kinetic term has to be provided. For this purpose, it is useful to 
restrict the arbitrary kets \fj,i >, i € M, to \fii >, i € i.e. to introduce the notion of regular graph 7^^, defined as a 
numerable set of equidistant points, whose separation is given by the parameter /zq, — {i ^ll = nfio,yn e Z}. 
The associated Hilbert space Ti-y^^ is separable. Because of the regular graph /xq, the eigenfunctions of must be of 



^ A Hilbert space is separable if and only if it admits a countable orthonormal basis. 
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the form e*™'*"'', m € 1i, which are Fourier modes, of period 2tt/^q. The inner product ([9]) is equivalent to the inner 
product on a circle with uniform measure, i.e.. 



< Hpmp) >Mo= ^ / Hp)m, (10) 



with p e (— TT/io, 7r//io), so that H-y^^^ = L^{S^,dp). Within this space, it is possible to construct an approximation 
for the "shift" operator, i.e. a regulated operator 

1 

PmoIM" >= 7^{\fJ-n+l > >)• (11) 

ZflQ 

More precisely, the polymer paradigm can be understood as the formal substitution 

—sininop), (12) 
Mo 

where the incremental ratio (jlip has been evaluated for exponentiated operators. The Hamiltonian operator iJ^o, 

which lives in Ti-^^.g, reads H^^ — + V{q), where the action of the new multiplication operator on wave 
functions in the momentum polarization is 



2 

Plo'^ip) = — [1 - cos(pAio)] , (13) 
Mo 

while the differential operator q is well defined. 

e. Continuum Limit. The physical Hilbert space of such theories can be constructed as the continuum limit 
of effective theories at different scales, and can be illustrated to be unitarily isomorphic to the ordinary one, 
ns^L^{R,dp). 

To this end, it is useful to remark that it is impossible to obtain Tig starting from a given graph 70 = 
{qk € R\qk = fcaojVfc g Z} by dividing each interval oq into 2" in new intervals of length a„ = ao/2", because Tis 
cannot be embedded into Hpoi- 

It is however possible to go the other way round and to look for a continuous wave function that is approximated 
by a wave function over a graph, in the limit of the graph becoming finer. In fact, if one defines a scale Cn, i.e., a 
decomposition of K, in terms of the union of closed-open intervals that have lattice points as end points and cover 
R without intersecting, one is then able to approximate continuous functions with functions that are constant on 
these intervals. As a result, at any given scale C„, the kinetic term of the Hamiltonian operator can be approximated 
as in (fT5|) . and effective theories at given scales are related by coarse-graining maps. In particular, it is necessary 
to regularize the Hamiltonian, treated as a quadratic form, as a self-adjoint operator at each scale by introducing a 
normalization factor in the inner product. The convergence of microscopically-corrected Hamiltonians is based on 
the convergence of energy levels and on the existence of completely normalized eigencoverctors compatible with the 
coarse-graining operation. 



IV. POLYMER TAUB UNIVERSE 



In this section, we will apply the polymer discretization technique to the description of the Taub model. In 
particular, we will specify the Hamiltonian ([8]) for the case of a discretized x space. As a result, the conjugate variable 
will not be implemented to operator directly, in the corresponding Hilbert space. Furthermore, the momentum space 
will be compactified, the compactification scale depending on the lattice characteristic length. 

The modifications to the Taub universe induced by the polymer representation will be investigated at both classical 
and quantum level. 



A. Classical Analysis 

First of all, let us clarify the physical meaning of our variables. The configuration variable x is related to the 
Universe anisotropy 7+ via the expression ©, for 61 = and ^ = (v^ + 3/4)/^/3w, as 
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FIG. 2: Semiclassical equations of motion for the Taub Universe: ordinary trajectory (blue line, a — 0) and polymer trajectories 
(green (cos a A — 1/2) and red {cos a A = 1/3) dashed lines). 



By this equation, a monotonic relation between the anisotropy of the Universe 7+ and our (classical) configuration 
variable x = Inv G [xq = ln(l/2),oo) appears, and, therefore, the variable x can be regarded as a measure of the 
model anisotropy. In particular, the isotropic shape of the Taub Universe (7+ = 0) comes out for a particular value 
of X, i.e. X = ln('\/3/2), and, in this case, we get the closed Friedmann- Robertson- Walker Universe. 
Let us now discuss the polymer dynamics of a Taub universe at classical level. By means of the substitution p^ . the 
Taub Hamiltonian reads 

H = — sin(a„p). (15) 
an 

From now on, we will take into account the discussion about the definition of a scale, and, for the sake of compact 
notation, we will drop the index n from a„. From a classical point of view, the equations of motion are 

X = {x, H} = cos(ap), (16a) 
p^{p,H}^0, (16b) 

where dot denotes differentiation with respect to the time variable t. The equations of motion are immediately solved 
as 



x(t) = cos(ap)T, (17a) 
p(t) = A, (17b) 

where A is a constant. 

As well understood, the system P7)) describes a free particle (Universe) bouncing against a wall. 
In the ordinary case, i.e. for a — 0, the model can be interpreted as a photon in the Lorentzian minisuperspace, and 
the classical trajectory in the (r — x) plane is its light-cone. More precisely, the incoming particle (r < 0) bounces on 
the wall (x = xq) and falls into the classical cosmological singularity (r —^ 00). 

Contrastingly, in the discretized case, i.e. for a 7^ 0, the one-parameter family of trajectories flattens, i.e. the angle 
between the incoming trajectory and the outgoing one is greater than tt/2 since p S (— 7r/a,7r/a) (see Fig. 2). As 
these trajectories diverge rather than converging, we expect the polymer quantum effects to be reduced with respect 
to the classical case, as we will verify below. 



B. Quantum Regime 

We now investigate the quantum behavior of the model. After analyzing the mathematical requirements of the 
polymer representation and their physical implications for the model, we apply the methods introduced above to the 
Taub Universe. In particular, we choose a discretized x space, and solve the corresponding eigenvalue problem in the 
p polarization. 

Even though the bulk of the discussion of the relation of the polymer representation at different scales is based on 
the properties of the Hamiltonian as a quadratic form, we can nevertheless apply this paradigm to the Taub model. 
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which is described by a Hnear Hamiltonian ([8]), after the well-known procedure, established in [21*1. In fact, squaring 
the Hamiltonian leads to squared eigenvalues without affecting the corresponding eigenfunctions. 

We are now ready to analyze the Schroedinger equation idr'i' — for the wave function ^' = r) corresponding 
to ([8]), where the configuation variable x is defined in the domain x € [xq = ln(l/2), oo). 

Considering the time evolution for the wave function as given by t) = e^^'''^ipk{p) and the results of [IH, we 

obtain the following eigenvalue problem 



2 

^ (1 — cos(ap)) — k 



Mp), (18) 



where, in the last step, the substitution (|13p has been taken into account. This eigenvalue problem is solved by 

= k\a) = 4 (1 - cos(ap)) < kl,, - 1 (19a) 

^k,a{p) = AS{p ~ pk,a) + B5{p + pk^a) (19b) 

ipk,ai^) = ^ [e:>q>iipk,ax) - exp(ipfe,a(2a;o - x))] : (19c) 

()19b[l is the momentum wave function, with A and B two arbitrary integration constant, and (|19cp is the coordinate 
wave function, where an integration constant has been eliminated by imposing suitable boundary conditions. Moreover, 
we have defined the modified dispersion relation 

Pk,a = - arccos ( 1 — j (20) 

from (|19ap . Furthermore, we stress that fc^ is bounded from above, as illustrated in (|19ap . but it is its square root, 
considered for its positive determination, which accounts for the time evolution of the wave function. 



V. TAUB WAVEPACKETS 



We will now gain insight onto the physical implications of the model by constructing suitable wave packets ^(x, r). 
In fact, analyzing the dynamics of such wave packets allows us to give a precise description of the evolution of the 
Taub model. Such an evolution will be preformed in both the polymer and Wheeler-DeWitt (WDW) approaches. 
More precisely, the latter will be considered the proper continuum limit of the polymer representation, as illustrated 
above. Wavepackets are a superposition of eigenfunctions (|19cl) . such as 

*(a;,r)=/ c^fcA(fc)^fc,a(a;)e-''=^ (21) 
Jo 

where A{k) is a Gaussian weighting function, i.e. A{k) — exp[— (fc — fco)^/2a-^]. 



A. WDW Dynamics 



To better understand the modifications induced on the ordinary dynamics by the polymer representation, we briefly 
summarize the WDW wave packet dynamics for the Taub model. In this case, the Hamiltonian is simply ([8]); the 
associated Schroedinger eigenvalue equation can be solved directly: the eigenfunctions in the position representation 
are just plane waves. This way, wave packets (1^ can be analytically calculated, with no upper limit for the energy 
k. The result is plotted in Fig. 3. As we can see from the picture, the wave packets follow the ordinary classical 
trajectories described in the previous Section. The probability amplitude to find the particle (Universe) is peaked 
around these trajectories. In this respect, no privileged regions arise, namely no dominant probability peaks appear 
in the (r — x) plane. As a matter of fact, the "incoming" Universe (r < 0) bounces at the potential wall (x — xq) 
and then falls towards the classical singularity (r oo). Therefore, as well-known, the WDW formalism is not able 
to shed light on the necessary quantum resolution of the classical cosmological singularity. As we will see below, this 
picture is slightly modified in the polymer representation. 
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FIG. 3: The WDW wave packet | ^(i, r) | for the Taub model, i.e. a = (fco = 0.1, a = 1). 




FIG. 4: The peaked polymer wave packet ] ^{x,t) \ for the Taub model, with koa = 1/2 (a = 50, fco = 0.01, cr = 0.0125). 



B. Polymer Dynamics 

We are now ready to analyze the modifications brought by the polymer representation in the quantized Taub 
Universe. Two cases can be distinguished, i.e the case k^a ~ 0(1), for which it is not possible to recover the 
ordinary representation of the momentum operator, and the case k^a <C 1, for which such a treatment is feasible. 
For /cqQ ~ C'(l), we get remarkable modifications of the wave packet evolution. From a probabilistic point of view, 
however, such modifications do not remove the cosmological singularity. The case k^a <C 1, contrastingly, can be 
considered as the semiclassical limit of the polymer approach. 

/. Peaked Weighting function. Let us now investigate the first case, k^a ^ 1, where the implementation of the 
polymer substitution (|12p does not lead to the ordinary Schroedinger dynamics. Furthermore, we stress that the 
choice of the value for the standard deviation a in the Gaussian weighting function can be relevant for detecting the 
effects of the polymer paradigm. 

In fact, if the weighting function is very sharply peaked around any value /cq, the resulting wave packet will be well- 
approximated by a purely monochromatic wave, for which a narrow neighborhood of k^ is selected. As a consequence, 
the ordinary dispersion relation is effectively reproduced by the deformed one, (j20p . In fact, narrowing the range 
of k is equivalent to expand the deformed Hamiltonian (jlSp around a given value of the momentum. This kind of 
behavior is explicitly illustrated in Fig. 4, where it is possible to appreciate a small interference phenomenon between 
the incoming (outgoing) wave and the wall. This feature can be interpreted as a relic of the polymer modifications of 
the Taub Universe dynamics, as it will be clearer in the next analysis. 

g. Spread Weighting function. On the basis of the previous analysis, the effects of the polymer substitution show 
up when broad wave packets are considered, i.e. when a large neighborhood of /cq is taken into account by the Gaussian 
weighting function. 
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FIG. 5: The spread polymer wave packet j ^{x,t) | for the Taub model, with koa — 1/2 (a — 50, ko — 0.01, a — 0.125). 



In this case, it is possible to appreciate all the modifications induced by the deformed Hamiltonian p^ . As a 
result, a strong interference phenomenon appears between the incoming (outgoing) wave and the wall. However, as 
a matter of fact, such an interference phenomenon is not able to localize the wave packet in a determined region of 
the configuration space. This way, the probability density to find the Universe far away the singularity is not peaked, 
i.e. the cosmological singularity of this model is not tamed by the polymer representation from a probabilistic point 
of view. Consequently, the incoming particle (Universe) is initially {t < 0) localized around the classical polymer 
trajectory (fl7|) . It then bounces against the wall {x = xo), where the wave packet spreads in the "outer" region, 
regains the classical polymer trajectory (t > 0) and eventually falls into the cosmological singularity {t — > oo). This 
way, we claim that the classical singularity is not solved by this quantization of the model. 

It is interesting to remark that the interference phenomenon occurs in the "outer" region of the configuration space, 
the (t — x) plane. These features are explained in Fig. 5. As we will discuss later on, such a behavior is complementary 
to that observed in the case of a generalized uncertainty principle. 



We end up our analysis by obtaining the correct semiclassical limit of the model. Within this framework, to 
obtain the proper continuum limit of the polymer representation, the value of ko is not arbitrary, but has to be 
chosen according to the request koa <^ 1. Since the range of the variable conjugated to the anisotropy variable is 
compactified, then ko has to be small with respect to the length of the interval*. 

As a result, differently from the other cases, the value of ko around which the wave packet is peaked is not arbitrary, 
but constrained by the characteristic scale a we are investigating. The ordinary WDW behavior is therefore recast, as 
plotted in Fig. 6. Even though taking ap <C 1 is enough to reproduce the ordinary Hamiltonian (as a general feature 
of the polymer representation because of relation p^ ). the fact that the correct semiclassical limit for the polymer 
quantum Taub Universe is obtained for a wave packet peaked at /cq ^ 1 is a non-trivial feature of the model. 



We can deeper understand the physical implications of this model by comparing it with other applications of the 
polymer representation in cosmology and with the implementation of a generalized uncertainty principle for the Taub 
Universe. In fact, in our model, the cosmological singularity is not probabilistically suppressed, as one could expect 
from other models. Let us now discuss the main differences from those models. 

h. Isotropic Polymer Cosmology. The fact that the cosmological singularity is not removed within this framework 
could look apparently in contrast with other models, such as [iJjUll- in the cosmological isotropic sector of General 
Relativity, i.e. the FRW models, the singularity is removed by loop quantum effects. In particular, the wave function 



We recall that the length of the integration interval L of II21I I is L oc 1/a, so that ko <JC L. 



c. 



Semiclassical Limit 



VI. COMPARISON WITH OTHER APPROACHES 



10 




FIG. 6: The semiclassical limit of the polymer wave packet | *I'(a;, r) | for the Taub model, with koa = 1/20 (a = 50, fco = 0.001, 
a = 0.01). 



of the universe exhibits a non-singular behavior at the classical sing ularity, and the Big Bang is replaced by a Big 
Bounce, when a free scalar field is taken as the relational time [23|. There are however at least two fundamental 
differences with respect to our model. 

• Within our scheme, the variable r, which describes the isotropic expansion of the Universe, is not discretized, 
but treated in the ordinary way. In fact, in the ADM reduction of the model, this variable emerges as the time 
coordinate, and cannot be discretized in a polymer frmework. More precisely, the phase space of this model is 
four-dimensional, but we naturally select a two-dimensional submanifold of it, i.e. the (x — px) plane, where 
we implement the polymer paradigm. In other words, we must discretize the anisotropy variable only, without 
modifying the volume (time) one. On the other hand, in the FRW case, the scale factor of the universe is 
directly quantized by the use of the polymer (loop) techniques. So far, the evolution itself of the wave packet of 
the universe is deeply modified by such an approach. 

• The solution of the equations of motion is radically different in the two cases. In fact, in our case, the variable 
p, conjugated to the anisotropy, is a constant of motion, and, from the Schroedinger equation, it describes also 
fc, the energy of the system. According to the polymer substitution p^ . it is always possible to choose a scale 
a for which the polymer effects are negligible during the whole dynamics, at classical level. On the other hand, 
the Hamiltonian constraint in the FRW case does not allow for a constant solution of the variable conjugated 
to the scale factor. For this reason, it is not possible to choose a scale, such that the polymer modifications are 
negligible throughout the whole evolution. 



i. Homogeneous Loop Cosmology. Also the Bianchi cosmological models have been analyzed in the framework 
of Loop Quantum Cosmology, according to the ADM reduction of the dynamics. The main difference between these 
works and our approach consists in the fact that in (2^ all the degrees of freedom are quantized by Loop techniques. 
In particular, also the time variable, i.e. the Universe volume, is treated at the same level as the others. In most 
cases, the time variable is defined by a phase space variable, i.e. it is an internal one. As a result, also the Bianchi 
Universes are singularity- free [25j . In this respect, our analysis is based on considering the time variable as an 
ordinary Heisenberg variable. 



]■ 



GUP Cosmology. The Taub universe, in this ADM reduction, has also been described within the framework of 



GUP [ll|. In that case, the conjugate variables x — px are quantized by means of a deformed Heisenberg algebra. As a 



result, the cosmological singularity is probabilistically suppressed, since the deformation parameter helps localize the 
wave function of the universe far away from it. This way, comparing the GUP approach and the polymer one allows 
us to infer that it i not always sufficient to " deform" the anisotropy variable to obtain significant modifications on the 
universe evolution. However, the polymer paradigm is a Weyl representation of the commutation relations, while, as 
explained in f26!|, a generalization of the commutation relations cannot by obtained by a canonical transformation of 
the Poisson brackets of the system. 

Moreover, it is possible to show how the effective framework of loop cosmological dynamics can be obtained by the 
opposite sign of the deformation term of the modified Heisenberg algebra [l^l . This feature is phenomenologically in 
agreement with our analysis. 
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VII. CONCLUDING REMARKS 



In this work, we have analyzed the polymer quantization of the Taub Universe. The Taub model admits a four- 
dimensional phase space, and its ADM reduction allows for an emerging time variable. So far, the energy variable and 
its conjugate momenta are treated canonically, while the anisotropy variable and its conjugate momenta are quantized 
according to the polymer paradigm. In particular, the anisotropy variable is assumed as discrete, while its conjugate 
momenta is replaced by its exponentiated version on a compactified space. 

This investigation has been developed at both classical and quantum levels. In the first case, trajectories are illustrated 
to flatten, with respect to the standard case. However, the most interesting result appears at the quantum level, when 
the evolution of wave packets is discussed. In fact, an interference phenomenon is illustrated to occur between the 
potential wall and the incoming particle (Universe), described as a localized wave packet. Nevertheless, the interference 
is not strong enough for the wave packet evolution to be localized. As a result, the corresponding outgoing particle 
(Universe) appears, whose evolution towards the cosmological singularity is not probabilistically avoided. 
The features of the polymer Taub Universe enhance the comparison with other approaches. On the one hand, the 
polymer quantization technique has been also applied to isotropic models. In this case, the choice of the scale factor 
as the polymer-discrete variable involves the singularity directly. This way, a non-singular quantum cosmology arises. 
On the other hand, the GUP approach to the Taub model leads to a singularity-free Universe. In particular, from an 
effective point of view, the consequences of the polymer scheme are complementary to those predicted by the GUP 
framework. 
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